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^ ; ABSTRACT 

The statistical mechanical investigation of violent relaxation (Lynden- 
Bell 1967) is extended to indistinguishable objects. It is found that, co- 
incidentally, the equilibrium distribution is the same as that obtained for 
classical objects. For massive neutrinos, the Tremaine & Gunn (1979) 
phase space bound is revisited and reinterpretated as the limit indicating 
the onset of degeneracy related to the coarse-grained phase space distri- 
bution. In the context of one of the currently most popular cosmological 
models, the Cold and Hot Dark Matter (CHDM) model (Primack et al. 
1995), the onset of degeneracy may be of importance in the core region 
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of clusters of galaxies. Degeneracy allows the neutrino HDM density to 
exceed the limit imposed by the Tremaine & Gunn (1979) bound while 
accounting for the phase space bound. 



Subject headings: dark matter — galaxies: clusters: general — methods: 
statistical 
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1. Introduction 

Self-gravitating astrophysical systems hav- 
ing initial mass and velocity distributions far 
from equilibrium are believed to evolve to- 
wards equilibrium by violent relaxation, a 
collective, collisionless process mediated by 
the fluctuating gravitational field. Compared 
to the time scale of the two-body relaxation 
process, the time scale of violent relaxation 
is short. The analytical investigation of vi- 
olent relaxation in terms of statistical me- 
chanics was first introduced by Lynden-Bell 
(1967). Formally, the equilibrium distribu- 
tion resembles the Fermi-Dirac distribution. 
For stellar and galactic systems it tends to- 
wards the Maxwell-Boltzmann distribution 
(Lynden-Bell 1967, Shu 1978). 

In this Letter we investigate violent relax- 
ation of indistinguishable objects. We show 
that, despite of the difference between the 
violent relaxation statistics of indistinguish- 
able and distinguishable objects, coinciden- 
tally, the equilibrium distributions of both 
cases are of the same form as expected by Shu 
(1987). When applied to massive neutrinos 
it turns out that the degeneracy implicit in 
the equilibrium distribution plays an impor- 
tant role leading to a reinterpretation of the 
Tremaine & Gunn (1979) phase space bound 
limit indicating the onset of degeneracy. 
In the context of one of the currently most 
popular cosmological models based on ini- 
tial conditions compatible with inflation the- 
ory (n = 1), the Cold and Hot Dark Mat- 
ter (CHDM) model (Primack et al. 1995), 
the onset of degeneracy allows the neutrino 
HDM density to exceed the limit imposed by 
the Tremaine & Gunn (1979) bound while ac- 
counting for the phase space bound. 



2. Violent Relaxation Statistics of In- 
distinguishable Objects 

Following Lynden-Bell (1967) we attempt 
to describe violent relaxation by the use of a 
maximum entropy principle. The one-particle 
phase space \x of the objects undergoing vio- 
lent relaxation is divided into a large num- 
ber of microcells the dimension of which is 
chosen so as to yield the mass density in the 
units mass per spatial and velocity volume, 
A 3 x and A 3 v, respectively. The collision- 
less nature of violent relaxation then imposes 
an additional constraint on the time evolu- 
tion of the /i-space. Initially not overlapping 
phase elements do never overlap. Therefore, 
a microscopic exclusion principle for phase el- 
ements in /x-space is established. The volume 
of the microcells is chosen in such a way that 
the microcells are occupied by at most one 
object each and that the probability of find- 
ing occupied adjacent microcells representing 
two-body encounters is small (Saslaw 1985). 
With m the particle mass, the mass per mi- 
crocell is either zero or m. Let r\ denote the 
corresponding mass density of an occupied 
microcell. The initial phase space distribu- 
tion before violent relaxation may then be de- 
scribed by the set of occupied microcells. 

At the macroscopic level, the microcells (or 
phase elements) are grouped into macrocells 
containing a large number v of microcells. 
Let N be the total number of occupied mi- 
crocells corresponding to N objects. Suppose 
that there is a microstate in which there are 
n a occupied microcells in the ath macrocell. 
Due to the collisionless interaction of the vi- 
olent relaxation process, the occupation state 
of each microcell does not change, i.e. there 
is no cohabitation. The indistinguishable na- 
ture of the objects under consideration leads 
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to a total number 



w a 



n„ 



W — n„ 



(1) 



of ways of assigning the n a occupied micro- 
cells to the v microcells of the ath macrocell 
as already realized by Lynden-Bell (1967) but 
not worked out further. Since the elements 
are indistinguishable, there is just one way of 
splitting the total of N elements into groups 
n a , i.e. there are no permutations. Thus, for 
indistinguishable elements, the total number 
of microstates is 



- A n a \(v- n a )\ 



(2) 



where, for distinguishable elements, the cor- 
responding number is (Lynden-Bell 1967) 



U a nJ a {y-n a )\ ' 



(3) 



The most probable state is found by the 
standard procedure of maximizing log W sub- 
ject to the constraints of constant total energy 
and constant mass. Since the expressions Wi 
and Wd differ only by a factor N\, the maxi- 
mization procedure performed in the contin- 
uum representation yields, coincidentally, for 
both expressions log Wi and log Wd the same 
coarse-grained phase space distribution 



f(v,x) 



n 



1 + exp[/5(e - fi)] 



(4) 



Here, e = v 2 /2 + $(x), with $(x) the gravita- 
tional potential normalized as $(x — > oo) = 
0. The Lagrange parameters /3 and \i are cho- 
sen according to the macroscopic constraints. 
Thus, the most probable coarse-grained phase 
space distribution of indistinguishable elements 
is the same as for distinguishable elements. 



Shu (1978) reexamined the statistical me- 
chanical discussion of violent relaxation in 
terms of particles and obtained an expression 
equivalent to (|J) determining the occupation 
number n a 



v 



n„ 



1 + exp[(3 p m(e a - /x p )] 



(5) 



of the ath macrocell where e a = v^/2 + ^(x a ), 
and the Lagrange parameters, (3 p and /x p , now 
refer to the particle description. For both the 
cases (d) and (H), with / <^ r/ and v a u, re- 
spectively, the phase space distribution tends 
towards the Maxwell-Boltzmann distribution. 
Despite of these agreements, the discussions 
by Shu (1978) and Lynden-Bell (1967) dif- 
fer what concerns degeneracy effects poten- 
tially present in ([|) and (|5|). While Lynden- 
Bell (1967) considers degeneracy to be even- 
tually important in central regions of galaxies, 
Shu (1978) concludes degeneracy to indicate 
the onset of two-body encounters. However, 
if two-body encounters are present, there is 
no microscopic exclusion principle any more. 
Thus the relevant phase space distribution is 
of Maxwell-Boltzmann type. In any case, for 
massive neutrinos, the discussion turns out to 
be irrelevant. 

3. The Case of Massive Neutrinos 

If the neutrino has a non-zero rest mass as 
predicted by some extensions of the Standard 
Model (Gelmini & Roulet, 1995), this may 
have important consequences for both the so- 
lar neutrino and the cosmological dark matter 
problem. Massive neutrinos with a rest mass 
of a few eV have been shown to account re- 
markably well for both problems (Primack et 
al. 1995). The discussion in this section is 
restricted to one neutrino flavor with g v = 1. 

We now apply the result of the foregoing 
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section to relict neutrinos of non-zero rest 
mass trapped in a gravitational potential, e.g. 
related to a cluster of galaxies. The initial 
phase space density of the unperturbed relict 
neutrino background is a relativistic Fermi 
distribution 



Up) 



m. 



(2nh) 3 eW^) + 1 



(6) 



with T v = 1.95 K (for z = and in the units 
mass per spatial and velocity volume d 3 x and 
d 3 v , respectively). For p = 0, the maximum 
value /„(0) = (l/2)ml/(2irh) 3 is obtained. 

If the neutrinos are massive, the phase 
space distribution / of the unperturbed relict 
neutrino background is changed due to gravi- 
tational interactions. Violent relaxation con- 
sidered here leads to a convolution of the 
phase space structure. Regions of initially dif- 
ferent phase space densities end up entwined 
together. However, the coarse-grained phase 
space density / cannot exceed the maximum 
of the initial fine-grained phase space den- 
sity f u since the collisionless nature of violent 
relaxation implies conservation of the fine- 
grained phase space density /. 

Since neutrinos are elementary particles 
and their unperturbed, initial phase space 
distribution is the relativistic Fermi distribu- 
tion (^j), the volume of the microcells is cho- 
sen to be the elementary phase space volume 
m 3 I \2irh) 3 . This choice guarantees the mi- 
crocells to be occupied by at most one ob- 
ject each. The corresponding phase space 
mass density of an occupied microcell is rj — 
ml/(2nh) 3 . Note, that the initial condition 
assumed in the foregoing section is perfectly 
fulfilled by the phase space distribution of the 
relict neutrinos. 

Suppose now the relict massive neutrinos 
to be subject to violent relaxation. According 



to the foregoing section, their coarse-grained 
phase space distribution / tends towards the 
Fermi type distribution (|]). First we note, 
that, in general, the relation max{/} < max{/} 
must hold. With the natural assumption 
P > 0, equation (01) implies 



($-//)> . 



(7) 



where for a given (3 equality corresponds to 
the highest possible phase space density. If 
no degeneracy effects are present, then / <C r\ 
and equation (H) requires 



P{Q-li) >0. 
In terms of the mass density 



(8) 



p„ = 4n V /Ve-^ 2/2+$ M^ (9) 



and the specified value of rj = mf, / \2nTi) 3 , in- 
equality (|8|) leads to 



4 3 
77V, (7 



V8^h 3 



(10) 



where a = (3 1//2 is the one-dimensional ve- 
locity dispersion. When solved for m v we get 




This expression is interpreted clS db limit on 
the type of the coarse-grained phase space 
density /. For given velocity dispersion a, 
neutrino density p v and a neutrino mass m v 
well above m^, / is the Maxwell-Boltzmann 
distribution. On the other hand, if m v is 
of order or even below m^, / differs from 
the Maxwell-Boltzmann distribution. In this 
case, degeneracy effects become important 
and the correct coarse-grained phase space 
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density / is the Fermi-type distribution (Q). 
We note, that the possible onset of two-body 
encounters preventing Fermi-like degeneracy 
in the case of stars or galaxies (Shu 1978) 
does not apply to neutrinos. The time scale 
of gravitational two-body encounters for neu- 
trinos is extremely long. Moreover, even 
when strong two-body encounters would oc- 
cur, the Fermionic nature of the neutrinos 
prevents single microcells to become cohab- 
itated. Thus, for neutrinos, a microscopic ex- 
clusion principle applies in any case. 

Inequality fllTD is formally equivalent to the 
well known Tremaine & Gunn (1979) bound 
on the rest mass of neutrinos which are gravi- 
tationally bound in galactic halos or galaxy 
clusters. However, the Tremaine & Gunn 
(1979) limit is derived from the assumption 
that / is the Maxwell-Boltzmann distribu- 
tion and the condition that the maximum of 
the coarse-grained phase space density / can- 
not exceed the maximum of the initial fine- 
grained phase space density / given by @. 
Therefore, the Tremaine & Gunn (1979) limit 
implicitly relies on the assumption of / to be 
of Maxwell-Boltzmann type. On the contrary, 
(ITTf) refers to the type of the phase space dis- 
tribution. Therefore we stress, that, from the 
point of view of ( |i~I"D , the limit implied should 
not be interpreted as limit on the neutrino 
mass. Instead, (|TT| ) states that the phase 
space distribution for m u near the Tremaine & 
Gunn (1979) bound differs significantly from 
the Maxwell-Boltzmann distribution. 

In terms of the gravitational potential $ 
and the neutrino mass m u it is possible to 
formulate several limits on the neutrino rest 
mass m v based on the phase space density 
/ defined by (|). We intend here to derive 
the most robust one. According to (|7|) the 
mass density of the totally degenerate, dens- 



est state is 



p v = J v dv 

o 



(12) 



where vi is the limiting velocity below which 
the maximum number of microcells is occu- 
pied (/ = 77/2) while above v\ all microcells 
are empty (/ = 0). Assuming only bound 
states to be relevant, i.e. e = v 2 /2 + $ < 0, 
the most robust limit on the neutrino mass is 
found for vi = 



m v > 



2$ I yielding 

fl2p u n 2 h 3 
\ |2$| 3 / 2 



1/4 



(13) 



This limit is robust in the sense that it holds 
for every coarse grained phase space distri- 
bution / representing a bound state includ- 
ing those eventually not described by as 
for example anisotropic phase space distri- 
butions suggested by Madsen (1991) or Ral- 
stone & Smith (1991). In order to compare 
these limits to the Tremaine & Gunn (1979) 
bound, the relation a 2 « $/3 holding for a 
Maxwellian velocity distribution related to a 
self-gravitating isothermal sphere is adopted. 
The Tremaine & Gunn (1979) bound on the 
neutrino rest mass then becomes 



m v > p 



1/4 



2\ 3/8 



(14) 



which is about a factor pa 1.18 higher than 
(P"B]). In terms of the neutrino density this 
leads to an augmentation of the density limit 
of about a factor of ~ 1.95. As a consequence, 
the onset of degeneracy allows the neutrino 
HDM density to exceed the limit imposed by 
the Tremaine & Gunn (1979) bound by a fac- 
tor of ~ 1.95 while accounting for the phase 
space bound. 
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In comparing the two limits, one should 
take into account that the Maxwellian ve- 
locity distribution related to ( |T4| ) has a tail 
of unbounded neutrinos. For the probably 
more physical truncated Maxwellian distribu- 
tion consisting of bounded neutrinos only, the 
difference between the limits would become 
bigger. On the other hand, the totally de- 
generate situation on which the limit flT3| ) is 
based is unlikely to be realized by violent re- 
laxation. 

4. Discussion and Conclusion 

This section considers violent relaxation of 
massive neutrinos in clusters of galaxies in 
the context of the currently popular Cold and 
Hot Dark Matter (CHDM) model. Exten- 
sive numerical simulations show these models 
to agree well with observations on cosmolog- 
ical scales (e.g. Klypin et al. 1993, Nolte- 
nius et al. 1994, Klypin & Rhee 1994). In 
what follows we refer to the two most promis- 
ing neutrino mass schemes due to Primack et 
al. (1995). The first scheme divides the neu- 
trino HDM in two neutrino species each with 
a mass of 2.4 eV while the second splits the 
HDM in three 1.6 eV neutrino species. The 
median value of the initial phase space dis- 
tribution f v of each family of the relict neu- 
trinos is only 0.06 (Madsen & Epstein 1984). 
As an approximation we therefore again take 
7] = m^/(27r/?,) 3 where m u is 2.4 eV or 1.6 eV, 
respectively, i.e. for our purposes the phase 
space elements occupied by two or three neu- 
trinos are neglected. 

The left panel of figure [I] presents the con- 
tour plot of density limits derived from ( |T3| ) 
and (|ToD for the 1.6 eV neutrino mass scheme. 
The right panel shows the density limits for 
the 2.4 eV mass scheme. The limits are plot- 
ted as upper limits to the total matter den- 



sities p consisting of HDM, CDM and bary- 
onic matter. They are functions of the grav- 
itational potential depth, indicated by the 
velocity dispersion of the corresponding self- 
gravitating isothermal sphere a p 



$/3, and 

the ratio of the neutrino HDM density p v to 
the total matter density p. The solid contours 
represent the most robust limit for bound 
states given by (|13|). Along these contours of 
specified total density p, the neutrino HDM is 
found in the totally degenerate, densest state. 
Augmentation of p is only possible if it is ac- 
companied by a change of the ratio p u /p or 
a. On the other hand, the more p u /p is di- 
minished while leaving p and a unchanged, 
the more the coarse-grained phase space / of 
the neutrino HDM becomes non-degenerate. 
The dashed contours are the density limit ob- 
tained from ([TDD . Well below this limit, / be- 
comes a Maxwellian. In addition, this density 
limit by ( |T4| ) corresponds to the Tremaine & 
Gunn (1979) bound. 

The parameter range of figure [I] includes 
typical values of a and p u /p expected for clus- 
ters of galaxies. The suppression of the HDM 
neutrino density p v compared to the CDM 
and baryonic matter density due to its lower 
density contrasts (e.g. Klypin et al. 1993, 
Noltenius et al. 1994) is reflected by values of 
p u /p below the cosmological value Q u = 0.2 
(Primack et al. 1995). For the 1.6 eV neu- 
trino mass scheme (left panel of figure [I]), 
one observes a broad parameter region yield- 
ing upper total density limits comparable to 
typical core densities of cluster of galaxies. 
Adopting the 2.4 eV neutrino mass scheme 
(right panel of figure [l]) shrinks the relevant 
parameter region. However, since the coarse- 
grained phase space distribution / of the neu- 
trino HDM is expected to be a Maxwellian 
only well below the limit indicated by the 
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dashed contours, it is found that for both 
CHDM models / may differ significantly form 
the Maxwell-type distribution. Accordingly, 
/ should be represented by the Fermi-type 
distribution (f|). Due to (|7|) it accounts for 
the phase space bound while it allows the neu- 
trino density to exceed the density limit im- 
posed by the Tremaine & Gunn (1979) bound 
(indicated by the dashed contour lines, too) 
by a factor of ~ 1.95. 

Violent relaxation eventually fades before 
the final state (§) is attained. Thus, @) is 
unlikely to be reached throughout the whole 
cluster, but it is reasonable to hold in the 
central region where violent relaxation occurs 
most violently. In order to determine the 
global properties of self-gravitating spheres 
based on coarse grained phase space distri- 
butions (f|), a careful investigation is needed. 
In particular, the onset of degeneracy even- 
tually related to core formation and not yet 
observed in simulations (Bryan et al. 1994) 
appears interesting. 

In summary, the statistical mechanical in- 
vestigation of violent relaxation has been ex- 
tended to indistinguishable objects. It is 
found that the equilibrium distribution is the 
same as that obtained for distinguishable ob- 
jects. However, in contrast to stellar and 
galactic systems, the onset of degeneracy is 
not prevented by two-body encounters in the 
case of massive neutrinos. Thus, the coarse- 
grained phase space density / is of Fermi- 
type and, to some extent, degeneracy may be 
present. For currently popular CHDM mod- 
els, the onset of degeneracy is shown to be rel- 
evant in the core region of clusters of galaxies. 

The authors thank Anatoly Klypin for valu- 
able criticism and discussion. 
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Fig. I. — Contour plot of upper limits to the total density p (in kg m~ 3 ) for the 1.6 eV (left 
panel) and the 2.4 eV (right panel) neutrino mass scheme. The limits are plotted as functions of 
the gravitational potential represented by the one dimensional velocity dispersion a according to 
a 2 Rs $/3 and the ratio p v j 'p. The solid contour lines represent the most robust density limit for 
bound states. Dashed contour lines indicate the density limit related to the type of the coarse- 
grained phase space distribution / of the neutrino HDM. Well below this limit / is a Maxwellian. 
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Fig. 1.- (left panel) 
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Fig. 1.- (right panel) 
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